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Abstract 

Two different solutions of the linearized Vlasov equation for finite systems, 
characterized by fixed and moving-surface boundary conditions, are discussed 
in a unified perspective. A condition determining the eigenfrequencies of 
collective nuclear oscillations, that can be obtained from the moving-surface 
solution, is studied for isoscalar vibrations of lowest multipolarity. Analytic 
expressions for the friction and mass parameters related to the low-enegy 
surface excitations are derived and their value is compared to values given by 
other models. Both similarities and differences are found with respect to the 
other approaches, however the close agreement obtained in many cases with 
one of the other models suggests that, in spite of some important differences, 
the two approaches are substantially equivalent. The formalism based on the 
Vlasov equation is more transparent since it leads to analytical expressions 
that can be a basis for further improvement of the model. 
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I. INTRODUCTION 



A collisionless kinetic equation with time-dependent mean field has been used long ago, 
first by Vlasov and then by Landau, to describe plasma oscillations in the high-frequency 
regime [Q. More recently Bertsch has proposed to use the same equation to study small- 
amplitude oscillations in nuclei [Q]. The approach of Vlasov and Landau assumes translation 
invariance of the system under study, extensions of their approach to finite systems have to 
face the nontrivial problem of which boundary conditions to impose on the fluctuation of 
the sigle-particle density distribution 6n{r,p,t). 

A solution of the linearized Vlasov equation suitable for describing giant resonances of 
various multipolarities in heavy spherical nuclei has been derived in The boundary 
conditions employed in [^] are inspired by the Steinwedel- Jensen (SJ) model of the giant 
dipole resonance (GDR) (compressible fluid within fixed boundary, see for example Ref. 
@], p. 558). In Ref. p an approach similar to that of Ref. P] has been followed, but with 
boundary conditions that are inspired by the Goldhaber- Teller (GT), rather than by the SJ 
model. In the GT model of the GDR the neutron and proton fluids are assumed to oscillate 
against each other without being compressed (see for example Ref. 0], p. 558). 

In order to adopt such a description, the authors of had to account for the reflection 
of nucleons on the moving nuclear suface (assumed to be sharp), and this was achieved by 
modifying the boundary conditions employed in 0. 

In this paper we first recall the approaches of Refs. and ||^ , giving a unified discussion 
of the different solutions derived in those papers. Then we consider the condition determining 
the eigenfrequencies of collective modes that has been derived in Ref. P]. Since in that paper 
this condition has been studied in detail only for monopole [L = 0) oscillations, we turn 
our attention to the lowest isoscalar modes with L > 0. We introduce a low- frequency 
approximation that allows us to obtain analytical expressions for the dynamic coefficients 
(friction and mass parameters) entering the equation of motion for the surface oscillations 
of a system of fermions. These coefficients and the associated eigenfrequencies of collective 
oscillations are then compared to the analogous expressions given by the liquid-drop model 
(LDM) of Ref. ^ and by a more closely related model discussed in Ref. 



II. ONE-DIMENSIONAL EXAMPLE 



In Ref. it has been shown that for a spherical system the solution of the linearized 
Vlasov equation can be reduced to the solution of a one-dimensional problem in the effective 
potential Uo{r) + 

In this section we briefly compare the approaches of and |^ by assuming a one- 
dimensional model, our aim is to make the connection between the two approaches more 
transparent. Thus we consider a one-dimensional system described by an equilibrium Hamil- 
tonian 

2 

/^o(x)=e=^ + f/o(x). (2.1) 

The (possibly self-consistent) equilibrium mean field U^ix) is assumed to have the shape 
of a potential well, so that the motion of particles is bounded within two classical turning 
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points Xi(2) defined by f/o(a;i(2)) = e. In the Hartree approximation the relation between the 
self-consistent field and the interaction between particles V{\ x — x' |) is 



Uq{x) = / dx'dp'V{\ x-x |)no(x',p'), (2.2) 



with no(x,p) the single-particle phase-space distribution at equilibrium. We assume that 
at time t = this system is acted upon by a weak external field of the type 5{t)Q[x) 
which induces a time-dependent fluctuation of the equilibrium single-particle distribution: 
riQ^x.p) — >• ?7.o(x,p) + 5n{x^p, t). As a consequence the equilibrium mean field also acquires a 
time-dependence: Uo{x) — ^ Uo{x) + 6U{x,t). It is convenient to take the Fourier transform 
with respect to time and to change variables from {x,p) to {x, e), so that 

6n{x, e, u) = 6n+{x, e, uj) + 6n_{x, e, uj), (2.3) 



with 6n±{x, e,Lu) = 6n{x,p = ±y2m[e — Uq{x)], lu). Then, according to the linearized 
Vlasov equation implies 

d5n+(x,e,uj) ., x ^ / \ 

A[x, e, uj)dn+[x, e, uj) = -d+(x, e, uj) 

ox 

d6n-(x,e,uj) x r / n n / n /o .x 

- -|- A[x, e, uj)dn_[x, e, uj) = B^[x, e, uj) (/.4j 



dx 



with 



and 



A(x,e,uj)= . (2.5) 



dno.JQ{x) dSU{x,uj) 
B±{x,e,uj) = (^)[^ + -^^]- (2-6) 



The solution of the system (|2.4| ) is not completely trivial because the fluctuation of the 
mean field 

6U{x,uj) = J dx'dp'V{\x-x' \)[6n+{x',e',uj) + 6n_{x',e',uj)] (2.7) 

couples the two equations. In the approach of P] the solution is obtained in two steps, 
first the mean field fluctuation in ( |2.6| ) is neglected, thus obtaining a zero-order solution 
5n^ that can be derived explicitly, then the solution with the full inhomogeneous term 
( p.6|) is expressed in implicit form through an integral equation similar to that given by the 
random-phase approximation (RPA) for the particle-hole Green function. 

Of course the zero-order solution 6n^ depends on the boundary conditions. The condi- 
tions employed in 0] were 

(5n+(xi(2),e,t^) = 5n_(xi(2),e,t^). (2.8) 
In this case the solution of the system ( ^.4[ ) (neglecting mean-field fluctuations) becomes 
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with B%x,e) = {^)m 



r{x) = f dy^=l== (2.10) 



e*^^ dx'B%x', e)e-*""(^') - /^^^^ dx' B\x' , e)e^'^"(^'') 

1 - e^' 



C(..) = v^^^>- ^ ^ ' ^ , (2.11) 



T = 2t{x2). (2.12) 
(We have simply re- written in more compact form the solution given in Sect. 4 of 0. Note 
that there is a misprint in Eq.(4.13b) of ||^, that should read: C2 = '^''~ilTxp{iuT) ^"^"^ )■ 

In Ref. it has been argued that the effects of moving boundary can be taken into 
account, at least in a linearized approximation, simply by changing the boundary conditions 
( p.8| ). In particular, since we are interested in the radial motion of nucleons in the effective 
potential f/o(r) + ^-j , it is sufficient to change the boundary condition at one turning point 
(say X2). Thus, instead of the boundary conditions (^.8]), in [^] the following boundary 
conditions have been used 

5n+(xi, e, oj) = 6n_{xi, e, u) 

5n+(x2, e, uj) = 6n_{x2, e, u) + (^)F(u;), (2.13) 

oe 

where F{uj) is the Fourier transform of an as yet undetermined function of time F{t). 
The solution analogous to ( p.9|) with these new boundary conditions is 

6h^ = 6n^ + 6n^s (2.14) 

with 

5ns{x,e,uj) = {—)F{uj) . . . 2.15 

oe zsm{ujT/2) 

We can obtain 5n'^{t) by taking the inverse Fourier transform of 5n^{uj) 

1 r+oo 

Sn\t) = — / duje~''^'6n\uj), (2.16) 

which can be determined by contour integration in the complex w-plane. Since we are 
interested in 6n^{t) for t > 0, the integration contour must be closed in the lower part of 
the complex plane. The first term on the r.h.s. of Eq. (|2.9|) is analytic in the whole plane, 
and the only contributions to the integral come from the poles of C{uj) that are at 

UJ = n— — ir] [n integer, rf — » ), (2-17) 

where 1 — e**^^ = 0. Equivalently we can thik of as a complex variable having a small 
positive imaginary part ir]. The moving-boundary solution is more complicated since the 
pole structure of 5n^g{uj) is not as simple as that of 5n^{ijj). Apart from the usual poles at 
UJ = n^, the quantity 6ng{u) can have additional poles due to F{lj), and these extra poles 
may contribute to the contour integral. The problem of determining the explicit form of 
F(t), and consequently the pole structure of F{uj), will be discussed in the next section. 
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III. THREE-DIMENSIONAL (SPHERICAL) CASE 



In this section we recall the approach followed in Ref. |^ in order to determine the pole 
structure of the additional term in the boundary condition ( |2.13| ). 

The approach of is inspired by the LDM the nuclear mean field is approximated 
by a square- well type potential, but, contrary to Ref. 0, in this model the external field 
is allowed to change also the equilibrium shape of the nucleus. The change of shape is 
parametrized as R{6, (p,t) = R + 6R{9, ip, t), with 

6R{e,<f,t) = J26RLM{t)YLM{0,f)- (3.1) 

LM 

In the LDM a surface deformation of this kind generates a fluctuation of the pressure at 
r = R that is given by (cf. Eq.(6A-57) of 0) 

6P{r,t) \r=R=Y.CL^-^^YLM{e,^). (3.2) 

LM ^ 

The restoring force parameter Cl used in [^] contains only the contribution of the surface 
energy, so 

Cl ^ {CL)surf = (yR\L - 1)(L + 2), (3.3) 
and a is the surface tension parameter that can be obtained from the mass formula: 

a ^ 1 MeV fm-l (3.4) 
In Ref. the link with kinetic theory was established by putting 

SP{r,t) \r=R= J dpprVr[Sn{r,p,t) - {^)SU{r,t)] \r=R . (3.5) 

This is a kind of self-consistency condition and its physical meaning is that the nuclear 
surface is required to behave as a free surface. The integral on the r.h.s. is the radial-radial 
element of the momentum flux tensor 6Ilrr{r,uj) for the nuclear quantum liquid p]. The 
radial velocity Vr{r) is 



(3.6) 

while the radial momentum is Pr = mvr- 

Combining Eqs.( |3.2|) and (|3.5| ), and taking the Fourier transform with respect to time 
gives 

aR''Y.iL-^)iL + 2)6RLM{uj)YLM{i)=SU„{r,uj) \r=R . (3.7) 

LM 

This equation can be used to determine the poles of the function F{lj) appearing in the three- 
dimensional generalization of Eq.( p.l3[) , that is, the eigenfrequencies of collective nuclear 
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vibrations. Multiplying both sides of Eg. ( |3.7| ) by y^j^^(r) and integrating over solid angle 
gives 

aR-\L - 1)(L + 2)6Rlm = OT.^f (r,a;) • (3.8) 

The r.h.s. is the multipole component of the integral in Eq.( p.5|) , it is given by a generaliza- 
tion of Eq.(B.6) of [§: 

{ lfe+(iA/JV,eAr.u-)+fc_(LMiV,eAr,^)]r;„(|,|) 

N=-L 

dno2L + l 



The quantities 5n±{LMN, e Xr,uj) satisfy the system of differential equations (|2.4|) with 



A{N,eXr,u) = ^ - ^4^^ (3.10) 
Vr{r) Vr[r)mr'^ 

replacing Eq.( |2.5|) and (there are some misprints in Eq.(5.12 of 0) 

B^iLMN, 6 A r, ^) = (^) [(|- ± ^^;^) (QlaKO + SU.Mir, ujWU^, f ) (3.11) 
instead of Eq.( p.6|) . Like in the one- dimensional example, the solution of the linearized 



Vlasov equation ( p.4|) can be studied in two steps. At first we can neglect, both in Eqs.( p.9|) 
and (|3.11|) , the mean- field fluctuation 6Ulm that couples the two equations, thus obtaining 
a zero-order approximation that may be a starting point for deriving the more complete 
solution that takes into account 6Ulm- In this paper we investigate only the zero-order 
approximation and show that even at this level some interesting results about low-energy 
isoscalar collective modes can be obtained. 

The zero-order solution analogous to Eq.(|2.9|) (fixed surface) is 



with 



J ri 

+ C{LMN,eX,u;)], (3.12) 



B^{LMN,eXr) = (^)[^ ± ^^;;^Q^m]>^^^( 2 2 )' (^.13) 



^(^) = / dy^—=—=r= , (3.14) 



A2 
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^(r) ^ r dy , ^ . ., , (3.15) 

C{LMN, eX,uj) = {e^[-^-^r]/2 / drBl{LMN, e A r)e-'[""(")-^^('-)] 

- r'dr5°(LMiV,eAr)e'[""(")-^^('')l}{l-e'["^-^^l}-\ (3.16) 

T = 2r(r2) T = 27(r2) . (3.17) 
The solution analogous to Eq. (|2.14|) (moving surface) instead, reads 
6fil{LMN, eXr,u) = 5n\{LMN, e\r,uj) 

+ (^)e^^l--(^)-^^(^)]F{LMN,eX,uj) -—-1—. (3.18) 



" V 2 ^ 

One more step is necessary in order to get the final result of Ref. 0: the functions 
F{LMN,eX,u!) must be related to the radial velocity of the surface (cf. Eq. (3.9) of [Q) 
through 

F{LMN,eX,Lu) = tu2prie X, R)6RLMiu;)YLNi^^)- (3.19) 



Replacing this expression into Eq.( p.l8D and recalling that 6n = Sn+ + (5n_, gives 6n^ = 
6n^ + Sn'g , with 

Sn's = 2(^) cos[urir) - Nj{r)] ujpr{tX, R)5Rlm{uj)Yln{11) ^^^^^r-NVy (3-20) 

Combining this equation with Eqs.(|^) and(|3.9|) and neglecting the mean-field fiuctuation 
^Ulm in (P-9|) gives the following expression for 6Rlm{^) in zero-order approximation: 



6R%,{uj) = (3.21) 
igL ^JdeJ dXXpr En Sn'jLMN, e X R, a;)y£^(f , f ) 
^(L^i^ - J de{^) J dXXpl En I YU^), f ) P cot(^^) ' 

As usual we assume that uj has a small positive imaginary part so that uj in the equation 
above has to be interpreted as uj + irj. The poles of 5R\j^.^{(jj) in the complex w-plane are 
determined by the vanishing of the denominator in Eq. (|3.21|) . This condition can be written 
in more compact form as 

Cl - Xl{uj) = (3.22) 

with Cl given by Eq.( p.3| ) and 

It is useful to recall that the functions lLAf(f , f ) vanish unless N has the same parity as L 
so that the sum over involves only terms with either odd or even A^. 

Eq.( |3.22| ) is the main result of Ref. ||^ and its solution determines the frequencies of 



collective nuclear excitations within that approach. 
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IV. COMPARISON WITH OTHER MODELS 



It is interesting to compare Eq.( |3.22| ) with the analogous condition given by the LDM. 



We consider the version that allows also for compression modes Q. Equation (6A-58) of 
Ref. can be written as 



Cl - /J% ] = (4-1) 



with Uc the velocity of sound. 

Equations ( [4.1|) and (|3.22| ) are coceptually similar because they both express a condition 
for the eigenfrequencies of collective nuclear modes. The main difference between them is 
that Eq. (4J.) has been obtained in a macroscopic hydrodynamic approach, while Eq.( p.22 ) 



has been derived within a microscopic kinetic-equation approach. It is thus reasonable to 
expect both similarities and differences between the eigenmodes determined by the two 
equations. 

For small u the term in square brackets in (|4.1| ) tends to 1 and Eq.( [4.1|) becomes 



Cl - l^'Dl = 0. (4.2) 

In this approximation the surface vibrations described by 6RLM{t) are pure harmonic oscil- 
lations and 5Rlm{'-^) has two simple poles at uo^ = ±-i/-§^ — if]- 



Dl 

In order to evaluate the function Xl{'^) "we take (^) = ~j^^^{^f — e), as appropriate 
for a gas of zero-temperature nucleons {ep is the Fermi energy, we use units such that % = 1), 



Even though Eq. (|3.22| ) is only a zero-order approximation, it is interesting since it allows 
us to obtain analytical expressions for the mass parameters Dl and for the other dynami- 
cal coefficients that determine the eigenfrequencies of collective modes. Unfortunately the 
functions Xiiyj) are not simple, only xl=o has been evaluated explicitly for any value of uj. 
The result derived in Ref. is, in terms of the dimensionless parameter s = 

°° 1 1 1 

XL=ois) = -2eFgoR'{l + 6 E[- + - - (4.3) 



n=l 



with 



W{s) = ]-s\n\^^\+t^sQi\s\-l), (4.4) 



^^o = the equihbrium density of nuclear matter and s„ = s/mr. 

We remark that, contrary to the hydrodynamic relation ( |4.1| ) that involves only real 
quantities, Eq. (p.22|) involves complex quantities and thus we can expect complex eigen- 



frequencies as solutions, in close analogy with the Landau damping phenomenon in homo- 
geneous systems. 

In this paper we want to extend the study of the functions xl{^), that was initiated in 
Ref. for the monopole case, to multipolarities L > 0. For this purpose we write Eq.( p.23D 
as 
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XLM = P|(o)xL=oM+Ar,.X: \ Yln{'^,^-) r/;vH, (4.5) 



with 



InIuj) = -uj 



/ (icosa cos a sin^ a [cot zjv(q;) + cot 2;_Ar(Q;)], (4.6) 
Jo 



2:7V « = ^ 75 + Na, 4.7 



and 



The quantities Pl{x) are Legendre polynomials of order L. It is useful to recall that 

^in„(f.f)P^p|(o)^o 

- ^ ■ 12 even L (4.9) 



^2^(f!)2 



To derive Eq. ([4.5| ) we have defined cos a = in Eq. (|3.23|) , the angle a has the simple 



geometrical interpretation shown in Fig.l. It follows immediately that Pr = PiT'Sina, T 



,T = 2a. 



/jv(c^) = 2cu / dxx^ — —7^ 77777 : r (4.11) 



2R sin a 

Vp 

The integrals 1^ can be complex because the integrands can have poles. In order to 
evaluate these integrals it is convenient to use the following trigonometric identity 

2 sinfzi + zo) , , 

cot zi + cot Z9 = z r 7 4. 10 

COS(Zi - Z2) - COS(Zi + Z2) 

and to change the integration variable to x = sin a, so that 

"1 - sin(2a:;s) 
/o cos[2a;(s + ii])] — cos(2A^arcsina;) 

The terms cos(2A^ arcsinx) = t/N are just polynomials inx^: yi = l — 2x^, y2 = 1— 8x^ + 8x'^, 
?/3 = 1 — 18x^ + 48x^ — 32x^, and so on. The integrals In can be easily evaluated for small 
s, thus we are led to study the following low-frequency expansion for the functions Xl{^)- 

Xl{^)^ AL + iu-iL + DLUj\ (4.12) 

This expansion is valid for s < 1, that is up to excitation energies of about MeV, 
however in the monopole case the upper limit of validity is vr times larger, since in that case 
the relevant expansion parameter is s/vr, thus we can expect to describe also the monopole 
compression mode within this approximation. 

The three coefficients A^, and given by the present approach can be compared with 
the same coefficients obtained in the hydrodynamic approximation and to those derived in 
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the more closely related approach of Ref. As we shall see, some of these coefficients vanish 
exactly, others (the mass parameters D^) may even diverge, in any case they convey non 
trivial information about the low-frequency behaviour of the functions XlI^^) cind SRlm{^)- 
The coefficients in Eg . (14.121 ) vanish for odd multipoles, while for even multipoles 

Al = -2epgoR'PliO). (4.13) 

The nonvanishing coefficients Al renormalize the restoring force parameter Cl in Eq.( p.22| ). 
The largest value of | Al \ occurs for L = when Ao = —2epQf)R^. For increasing (even) 
L the coefficients Al tend to zero as l/L, while in Eq. ( |3.3| ) increases as L^, so this 
renormalization becomes negligible in the limit of large L. However for L = A/, is much 
larger than Cj^, in agreement with the well known LDM result that the monopole excitation 
is not a surface, but a compression mode (in this case the surface tension plays a negligible 
role). For L = 2 the value of Al is comparable to that of Cl and this leads to a reduction 
of the role played by the surface tension in this excitation mode too. 

The friction coefficients 7^ are most easily evaluated directly from Eq.( |3.23D by using 
the pole expansion 

+00 

cotz= y . (4.14) 

It is clear also from Eq. (|4.4|) that for L = the friction coefficient vanishes, while for L > 
we find 

(47r)^ ^1 TT TT TT 

7l = 7wf ^ 2L + l ^ iV ' ^^^''2' 2'' ^ cosa„jv sin^ a„ivO(- - ^un) (4.15) 

with anN = -^TT. The angles UnN appearing in this equation are related to the nucleon 
trajectory in the unperturbed mean field in the way shown in Fig.l . It is interesting to note 
that when the coefficient Jl is nonvanishing, it gets a contribution only from nucleons moving 
along closed classical trajectories. Because of the vanishing of the coefficients Yln{^, f ) 
when the parity of differs from that of L, for even multipoles the sum over in Eq.( [4.15D 
effectively starts from N = 2 and involves only even values of A^. 

For ease of comparison with Ref. we have expressed our values of 7^ in terms of 
the wall formula friction coefficient 7^;/ = ^QqPfR'^- We find that Eq.( |4.15| ) reproduces 
very well the values of 7^ reported in Ref. 0. For all value reported there (up to L = 12) 
we find perfect agreement for odd multipoles, while for even multipoles our values of 7^ 
are systematically about 3 — 4% smaller than the values of 0. This discrepancy tends to 
decrease with increasing L, and is only about 1% for L = 12. 

From Eqs.( |4.5| ) and ( [4.6| ) we can obtain an explicit expression for the coefficients -Dl in 
the expansion (^l2D , it reads 



1 AtT 71 71 

Dl = -mgoR' Pl{0) + ^^^^oR'^^J2 I Yln{-, -) T 4 (4-16) 

with 
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Note that possible terms that behave as 1/s resuhing from the evaluation of I'^ should be 
omitted from Dl since they have already been included in the coefficients 7l and that the 
sign of the imaginary part of depends on the sign of s. 

The first term on the r.h.s of Eg. ( |4.16D is present only for even multipolarities and is 



the only term for L = 0, the remaining sum over involves only even or odd values of N ^ 
depending on the parity of L. 



The first few integrals are: 



= \. (4.18) 



11 7r 
IL = — lim In I s I +-(2 In 2 - 1) + i — sign s (4.19) 



1^3/.^ , 1 + ^3/2, 7rV3, 1 , 

I' = — 4V3+n -^)+i——\im- 4.20 

^ 64 3 ^ 1-^3/2^ 16 2 ^-os ^ ^ 

For all even multipoles with L > the real parts of the integrals are divergent because 
there is a pole exactly at the upper integration limit. This divergence agrees with that 
found in Ref. In the present formalism it is possible to pinpoint exactly the origin of 
this divergence: it is generated by nucleons moving on orbits along the diameter. A peculiar 
feature of the present approach is that the coefficients Dl can be complex. In this aspect 
our results differ radically from those of Refs. and 0. Our coefficients Dl are real for 
L = 0, 1, but for L = 2 the real part is divergent and the imaginary part is finite. Note that 
this does not imply that the function Xl=2 has a pathological behaviuor for u ^ 0, since 
the divergence is only logarithmic. 

The integral I'^ is also complex, however its imaginary part behaves like 1/s, the corre- 
sponding term in the expansion (|4.12|) is linear in u and it has already been included in the 



coefficient 7^. Thus for example for L = 3 only the real part of I'^ should be included in 

In Table 1 we compare the values of the dynamic parameters given by the present ap- 
proach with those of Refs. for the first few multipoles. 

As already found in Ref. the present approach gives the following value for the energy 
of breathing modes in nuclei: 



hujL=o ~ Hh^jT^ ~ 96 A-^/^ MeV. (4.21) 
V -Do 

This value can be compared with that of 65 A~^^^ MeV reported in Ref. for classical 
hydrodynamics (both numbers have been obtained by assuming the parametrization R = 
1.2 A^^^ fm, for this reason our value differs from the value 103 A~^^^ MeV of Ref. |^). 
We recall that the result ( [4.21|) is based on a zero-order approximation in which mean-field 
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fluctuations in the bulk have been neglected, taking into account these effects leads to fTO 
hujQ ~ 82 A~^^^ MeV, in reasonable agreement with experiment [|Tl|] . 

The L = 1 case is the only one in which all models of Table 1 give the same results. This 
simply means that none of them has problems in describing centre-of-mass motion. 

For L = 2 the eigenfrequency condition ( |3.22| ) cannot have solutions corresponding to 
real, small, u because the coefficient -Dl=2 is complex. Thus we have to look for possible 
solutions in the complex uj or s plane. For this puspose we write Eq.( ^.23| ) as 

Xl=2('S) = 3Aos{- / dxx^ cot(sx) + - dxx^[cot{sx — 2a) + cot(sx + 2a)]}, (4.22) 
4 JO 8 JO 

with a = arcsinx. For small s, cot(sx) ~ l/sx — sx/3 + O(s^), and the first integral is easily 
evaluated. Similarly, the second integral, using the identity (|4.1CI| ), becomes 



/ dxx^[cot{sx — 2a) + cot{sx + 2a)]} 
Jo 



1 x'' 



4s / dx-—: 7^—-^ + 0{s 



k 8(x2 - 1) + 2s2 

-2^2^''^^ - S 
Thus for complex s and |s| < 1 we have 



-o[oMl--)-l] +0{s% (4.23) 



If we use this approximation for X2(s) the eigenfrequency condition (|3.22|) has no so- 
lution with \s\ < 1 in the complex s-plane, thus the present model, at least in its simple- 
minded zero-order approximation, does not seem to be adequate for describing low-frequency 
isoscalar quadrupole oscillations in nuclei. Calculations of the isoscalar quadrupole response 
performed in Ref. |ll2|] within the same model used here (for real u) suggest the existence 
of a collective mode in a region where |s| ~ 1 (cf. Fig. 2b of [0). Taking into account the 
attractive interaction inside nuclei (that is the term 5Ulm in Eqs. (|3.9| ) and (|3.11|) ) could 
bring this collective mode into the region |s| < 1, however this would require a more numer- 
ical approach and would mean loosing the insight given by the analytic expressions of the 
zero-order approximation. 

For L = 3 the eigenfrequency condition ( |3.22| ) in the limit of small uj becomes 

ia;73 + D^uj^ = C3. (4.25) 

Using the values of parameters reported in Table 1, we find that low-energy octupole oscilla- 
tions are overdamped. This result might offer a qualitative explanation for the background 
observed in inelastic proton scattering at small angle, that has been previously studied by 



assuming a semi-infinite model |13|, |14 



V. SUMMARY AND CONCLUSIONS 

Different solutions of the linearized Vlasov equation for finite systems can be obtained, 
depending on the boundary conditions imposed on the fluctuations of the single-particle 
phase-space density. 
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We have discussed two different solutions (with fixed and moving surface) that can 
be useful in studying collective nuclear excitations. In particular we have analyzed the 
moving-surface solution in the low-frequency limit for the first few isoscalar nuclear modes 
(monopole, dipole, quadrupole and octupole) and have derived interesting analytical ex- 
pressions for the friction and mass parameters entering the equation of motion for the os- 
cillations of the nuclear surface. According to our semiclassical result, in the low-frequency 
limit, nuclear dissipation is due only to nucleons moving along closed classical trajectories. 
Numerically our friction parameters are in rather good agreement with those derived in Ref. 
|[^. For the mass parameters our result agree only partially with those of ||^, like in that 
case we find a divergence of the mass parameters for even multipolarity, but our parameters 
can be complex, with a finite imaginary part. For odd multipolarities (dipole and octupole) 
our mass parameters do agree with those of Ref. 0]. Thus, even though the formalism is 
quite different, we conclude that the approach of Ref. is substantially equivalent to that 
of Ref. [0 . On the contrary we find several differences with respect to the classical hydrody- 
namic results of Ref. p . These differences, the most imporant being probably the Landau 
damping of surface oscillations, can be ascribed to the fact that the properties of a fluid of 
fermions are, even in a semiclassical description, different from those of a classical fluid. 
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FIGURES 




FIG. 1. Relation of angle a to nucleon trajectory in equilibrium mean field 
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TABLES 



TABLE I. Dynamic parameters for lowest multipolarities. 
Columns a: present results, columns b: Ref.[7], columns c: model of Ref.[6] 
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